OFFICE  OF  NAVAL  RESEARCH 

GRANT:  N00014-90-J-1230 

R&T  CODE  4133015 

ROBERT  J.  NOWAK 

Technical  Report  No.  63 

Electron  Transfer  Kinetics  of  Self-assembled  Ferrocen  C(12)  Alkanethiol  Monolayers  on  Gold 

Electrodes  from  125K  to  175K 

by 


John  N.  Richardson,  Gary  K.  Rowe,  Michael  T.  Carter,  Leonard  M.  Tender,  Larry  S.  Curtin, 

Stephen  R.  Peak  and  Royce  W.  Murray 


Prepared  for  Publication 
in  the 

Electrochimica  Acta 

University  of  North  Carolina  at  Chapel  Hill 
Department  of  Chemistry 
Chapel  Hill,  NC 

March  1,  1995 

Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 

United  States  Government 

This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is  unlimited. 


DTIC  wi/AiJxT  mnpBOTED  1 

1 9950328  040 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0738 

Puanc  reDomng  burden  for  :nis  collection  of  inrormation  is  estimaiec  to  average  ’  hour  per  response,  including  the  time  for  reviewing  instruaions.  searcning  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  anv  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burcen.  to  Washington  Headauarxers  Services.  Directorate  for  information  Operations  and  Reports  1215  iefferson 
Davis  Highway.  Suite  1204,  Arlington.  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduaion  Project  (0704-0 188).  Washington.  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

March  23,  1995  Interim 

4.  TITLE  AND  SUBTITLE 

Electron  Transfer  Kinetics  of  Self-Assembled  Ferrocene  C(12) 
Alkanethiol  Monolayers  on  Gold  Electrodes  from  125K  to  175K 

5.  FUNDING  NUMBERS 

N00014-90-J-1230 

6.  AUTHOR(S) 

lohn  N.  Richardson,  Gary  K.  Rowe,  Michael  T.  Carter,  Leonard  M. 
Fender,  Larry  S.  Curtin,  Stephen  R.  Peak  and  Royce  W.  Murray 

7.  PE.RFCRMING  ORGANIZATION  NAME(S)  AND  ADORE5S(ES) 

The  University  of  North  Carolina  at  Chapel  Hill 

Department  of  Chemistry 

Chapel  Hill,  N.C.  27599-3290 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

63 

9.  SrCNSORl.NG/ MCNITCaiNG  AGENCY  NAME^S)  AND  ADDRESSEES) 

Office  of  Naval  Research 

Department  of  the  Navy 

Arlington,  VA  22217 

10.  SPCNSCRING/MONITCRING 

AGE.NCY  REPORT  NUMBER 

j11.  SUPPLEMEI'-iTARY  NCTES  i 

I  .  I 

1:3.  C.-i  -711  i~,->  rS.vlEN  ■ 

>  20~  0 i b  1  ,Ri  i u,  .  ;  0 1'f  C C  U"  d  1 

1 

Approved  for  Public  Release,  Distribution  Unlimited 


13.  ABSTRACT  [Maximum  2C0  wcrazi 


The  electron  transfer  kinetics  of  monolayers  of  CpFeCpC07(CH2)i2SH  and  CH3(CH2)iiSH  co- 
chemisOTbed  on  gpld  electrodes  have  been  measured  in  2:1  (v:v)  chloroethane/butyronitrile  solvent  at 
temperatures  ranging  from  125K  to  175K  with  potential  steps  and  cyclic  voltammetry.  Rate 


constants,  k‘^ 


measured  using  cyclic  voltammetry  range  from  3x10^  to  10’'  s'^ 


over  these 


temperatures;  an  activation  plot  of  log[k7kbT'^]  vs.  1/T  gave  X  =  0.89  eV  for  the  reorganization 
energy  for  the  Cp2Fe'^™  reaction  and  a  pre-factor  of  6.5x10^  eV'^s  ^  The  rate  constants  were  obtained 
by  comparison  of  experimental  AEpeak  values  with  those  of  cyclic  voltammograms  digitally  simulated 
with  Marcus-DOS  theory.  A  value  of  X  =  0.76  eV  is  predicted  from  dielectric  continuum  theory. 
The  CpFeCpC02(CH2)i2SH  monolayers  are  kinetically  inhomogeneous  (disperse)  as  shown  by  non¬ 
linear  ln[i]  vs.  time  plots  in  potential  step  experiments.  The  kinetic  dispersity  has  adverse  effects  on 
determinations  of  X  values  from  cyclic  voltammetric  waveshapes  and  from  plots  of  potential  step- 
derived  rate  constants,  k^pp  ,,  against  overpotential  {rj),  producing  depressed  values  for  X. 


14.  SUBJECT  TERMS 

Ferrocene,  Electron-transfer  kinetics.  Monolayers 


17.  SECURITY  CLASSiFICATlON 
OF  REPORT 

Unclassified 


18.  SECURITY  CLASSIFiCATlON 
OF  THIS  PAGE 

Unclassified 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 


15.  NUMBER  OF  PAGES 

32 


16.  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 

Unlimited 


MSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  $td.  Z39-'8 


t 


October  12,  1994 

ELECTRON  TRANSFER  KINETICS  OF  SELF-ASSEMBLED  FERROCENE 
(C12)ALKANETHIOL  MONOLAYERS  ON  GOLD  ELECTRODES  FROM  125K  to  175K 
John  N.  Richardson,  Gary  K.  Rowe,  Michael  T.  Carter,  Leonard  M.  Tender^,  Larry  S. 
Curtin^,  Stephen  R.  Peck^,  and  Royce  W,  Murray 

Kenan  Laboratories  of  Chemistry,  University  of  North  Carolina  at  Chapel  Hill,  Chapel  Hih, 
NC  27599-3290 

ABSTRACT 

The  electron  transfer  kinetics  of  monolayers  of  CpFeCpC02(CH2)i2SH  and 
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of  experimental  AEpeak  values  with  those  of  cyclic  voltammograms  digitally  simulated  with 
Marcus-DOS  theory.  A  value  of  A,  =  0.76  eV  is  predicted  from  dielectric  continuum  theory. 
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by  non-linear  ln[i]  vs.  time  plots  in  potential  step  experiments.  The  kinetic  dispersity  has 
adverse  effects  on  determinations  of  X  values  from  cyclic  voltammetric  waveshapes  and  from 
plots  of  potential  step-derived  rate  constants,  Kpp,r,^  against  overpotential  (77),  producing 
depressed  values  for  X. 
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ELECTRON  TRANSFER  KINETICS  OF  SELF-ASSEMBLED  FERROCENE 
(C12)ALKANETHIOL  MONOLAYERS  ON  GOLD  ELECTRODES  FROM  125K  to  175K 

Molecular  monolayers  organized  by  hydrophobic  forces  are  long-standing  objects  of 
research  in  chemistry  and  molecular  biology.  In  recent  years,  the  combining^  of  hydrophobic 
forces  and  surface  chemisorption  has  opened  new  avenues  for  fabrication  and  study  of 
molecular  films.  In  particular,  chemisorption  and  self-assembly  of  alkanethiols  and  their 
derivatives  on  metal  surfaces,  notably  gold,  is  undergoing  detailed  scrutiny  in  investigations 
of  chemically  modified  electrode  interfaces^  and  of  the  dynamics  of  electron  transfers 
between  the  electrode  and  reactants  positioned  on  the  solution  side  of  an  organized  alkane 
monolayer  barrier.^  There  have  been  important  corollary  studies  on  mixed  monolayers'*  and 
monolayer  structure^,  and  related  investigations  of  chemically  modified  metal  surface 
wetting®,  corrosion^,  and  patterning®.  Measurements  of  electrode  kinetics  at  self-assembled 
monolayers®  promise  to  provide  important  new  insights  into  how  the  detailed  molecular 
environment  of  an  electrode  reactant  affects  its  reactivity®'’,  and  into  the  chemistry  of 
electron  tunneling®.  Control  of  the  distance  between  electrode  reactants  and  electrode  was 
a  long-standing  goal  in  research  on  chemically  modified  electrodes®",  and  has  become  a 
recent  target  in  design  of  interfaces  for  study  of  electron  transfer  dynamics  at 
superconducting  electrodes.®^ 

This  paper  will  describe  measurements  of  low  temperature  electron  transfer  kinetics 
in  self-assembled  mixed  monolayers  of  ferrocene  ester  dodecanethiol, 
CpFeCpC02(CH2)i2SH,  and  dodecanethiol,  CH3(CH2)iiSH.  The  Cp2Fe^/"  electrode 
reaction  is  studied  in  a  binary  cryoelectrochemical  solvent®®,  2:1  (vrv) 
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chloroethane:butyroiiitrile  (Eta/PrCN)  over  a  temperature  range  of  125K  to  175K, 
measuring  its  standard  electron  transfer  rate  constant  k"  using  cyclic  voltammetric  and 
potential  step  experiments.  We  have  previously  shown^  that  CpFeCpC02(CH2)i2SH  and 
shorter  chain  CpFeCp(CH2)3SH  monolayers  are  exceptionally  stable  in  this  low  temperature 
medimn.  Experiments  on  the  electrode  kinetics  of  CpFeCp(CH2)gSH  and  of 
CpFeCpC02(CH2)8SH  monolayers  on  Au  and  Ag  electrodes  are  described  elsewhere.^ 
Central  questions  in  the  present  study  have  been  the  suitability  of  various  approaches  to 
measurement  of  k°  and  of  the  reorganization  energy  X,  and  whether  these  parameters  exhibit 
any  unusual  properties  as  a  result  of  the  lowered  temperatures  or  of  the  organic  solvent 
medium.  These  questions  are  inter-connected  because,  in  contrast  to  those  described  by 
Chidsey  et  al^,  the  ferrocene  ester  alkanethiol  monolayers  studied  here  are  not  kinetically 
homogeneous,  probably  in  part  due  to  the  organic  solvent  environment.  The  kinetic 
inhomogeneity  has  adverse  consequences  on  the  suitability  of  different  measurement 
approaches  to  X. 

EXPERIMENTAL 

Chemicals.  Ferrocene  ester  dodecanethiol,  CpFeCpC02(CH2)i2SH  (Cp  =  cyclopentadienyl) 
was  prepared  by  a  literature  method.*^  Dodecanethiol  (CH3(CH2)iiSH,  Aldrich  97%), 
butyronitrile  (Aldrich,  99+  %),  absolute  ethanol  (AAPER  Alcohol  and  Chemical  Co.),  and 
tetra-n-butylammonium  hexafluorophosphate  supporting  electrolyte  (BU4NPF6,  Fluka,  Puriss) 
were  used  as  received.  Chloroethane  (Linde)  was  condensed  and  stored  in  sealed  vacuum 
transfer  pipets  at  room  temperature.  Water  was  purified  ( >  18  Mohm/cm^)  with  a 
Bamstead  NANO  pure  system. 
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Electrode  Fabrication  and  Preparation.  Working  and  reference  electrodes  were  0.5  mm  and 
2.0  mm  diameter  gold  rods  (Johnson  Matthey  Electronics,  99.999  %),  respectively, 
encapsulated  in  and  exposed  at  the  end  of  an  epoxy  cylinder  (Shell  Epon  828,  m- 
phenylenediamine  curing  agent,  cured  overnight  at  70°  C).  The  metal  disks  were  polished 
before  each  experiment  with  1/im  alumina  (Buehler)  and  1  jum  diamond  paste  (Buehler) 
with  extensive  rinsing  with  water  and  sonication.  Prior  to  immersion  in  the  thiol  coating 
solution,  the  gold  surfaces  were  etched  in  dilute  aqua  regia  for  five  minutes^^  and  thoroughly 
rinsed  with  water  and  ethanol. 

Chemisorbed  Monolayer  Preparation.  Coating  solutions  consisted  of  a  1:3  (M:M)  ratio  of 
CpFeCpC02(CH2)i2SH  and  CH3(CH2)iiSH  in  ethanol  at  ImM  total  thiol  concentration. 

The  CH3(CH2)iiSH  diluent  thiol  is  employed  to  minimize  lateral  interactions  between  the 
ferrocene  sites'*^.  The  etched  electrodes  were  exposed  to  this  solution  for  >24  hours, 
thoroughly  rinsed  with  ethanol,  and  dried.  The  Au  pseudo-reference  electrode  was  scraped 
to  expose  a  fresh  Au  surface.  Annealing-incubation  of  chemisorbed  films  in  fresh  diluent 
solution  was  not  employed.**® 

Electrochemical  Measurements.  The  cylindrical  electrode  assembly  was  fitted  into  a  slotted 
stainless  steel  sleeve  whose  Pt  bottom  served  as  an  auxiliary  electrode  spaced  ca.  0.5  mm 
from  the  working/reference  electrodes.  The  electrodes  were  placed  in  an  aluminum 
container  containing  the  desired  quantity  of  butyronitrile  and  BU4NPF6  (0.075M  final 
concentration),  pre-chilled  to  dry  ice  temperature  to  avoid  thiol  desorption.  The  cell  was 
sealed,  immersed  in  a  dry  ice/acetone  slush  bath,  evacuated,  and  chloroethane  was  vacuum 
transferred  to  the  final  solution  volume.  The  cell  was  bolted  to  the  cold  finger  of  a  Janis 
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helium  refrigerator  cryostat  interfaced  to  a  Lakeshore  Cryogenics  320  Autotuning 
Temperature  Controller,  and  brought  to  the  desired  temperature,  which  required  ca.  4  hour. 
Further  temperature  changes  of  <10°  required  ca.  1  hour  for  equilibration. 

Cyclic  voltammetric  and  potential  step  experiments  were  executed  with  instruments  of 
local  design,  employing  A/D  conversion  and  digital  data  storage,  low-pass  filtering  of  analog 
signals  during  slow  timescale,  and  positive  feedback  iR  compensation  to  counter  the  low 
solution  ionic  conductivities.  The  latter  step  was  particularly  important  given  that,  for 
example,  without  feedback,  the  potentiostatically  uncompensated  resistance  between  working 
and  reference  electrodes  at  125K  was  as  large  as  550kn.  (Uncompensated  solution 
resistances  were  measured  by  AC  impedance  with  a  Schlumberger-Solartron  Instruments 
Model  1255  frequency  response  analyzer  and  Model  1286  potentiostat.)  Positive  feedback 
in  cyclic  voltammetry  and  potential  step  experiments  was  implemented  by  incrementally 
increasing  the  proportionality  constant  between  cell  current  and  feedback  potential  added  to 
the  applied  electrode  potential  until  the  circuit  became  slightly  unstable  (current  oscillation), 
then  backing  off  to  restore  circuit  stability.  We  have  found  that  such  instrumental 
correction  for  resistance  effects  is  well-suited  to  low  ternperamre  work  in  which  the 
measured  currents  are  small  and  experimental  time  scales  long. 

Potential  step  experiments  were  initiated  from  the  formal  potential  (E°’,  determined 
by  cyclic  voltammetry)  of  the  Cp2Fe^'^°  couple;  the  applied  overpotentials  varied  from  ±30 
mV  to  ±500  mV.  Corrections  for  background  currents,  which  were  minor,  were  obtained 
from  potential  steps  in  the  double-layer  potential  region  negative  of  E° .  Coverages  of 
electroactive  ferrocene  (r)  were  measured  from  the  charge  under  slow  potential  scan 
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(lOmV/s)  cyclic  voltammograms. 

RESULTS  AND  DISCUSSION 

Measurements  of  Heterogeneous  Electron  Transfer  Rate  Constants  With  Cyclic 
Voltammetry.  Typical  examples  of  low  temperature  cyclic  voltammograms  of  mixed 
CpFeCpC02(CH2)i2SH/CH3(CH2)iiSH  monolayers  are  shown  in  Figure  1.  The 
voltammograms  are  well-defined,  show  symmetrical  oxidation/reduction  branches,  and  are 
extremely  stable,  remaining  unchanged  over  hours  of  potential  cycling  at  these  temperatures. 
The  quantities  of  ferrocene  sites  present  in  the  monolayers  are  governed  by  the 
ferrocene/diluent  ratio  in  the  coating  solutions;  typical  ferrocene  coverages  (Figure  1)  are 
ca.  1/4  to  1/3  of  a  full  monolayer.'*^’^^ 

Figure  1  shows  that  the  separation  AEp^AK  between  oxidation  and  reduction  peak 
potentials  in  CpFeCpC02(CH2)i2SH  monolayer  voltammetry  increases  as  the  electron 
transfer  rate  is  thermally  quenched  by  lowered  temperatures.  AEpE^K  ^so  increases  at 
increased  potential  sweep  rate  as  shown  in  Table  I. 

From  theory  by  Laviron^'*  based  on  the  Butler-Volmer  formalism  for  electrode 
kinetics^,  AEp^AK  is  governed  by  the  electron  transfer  rate  constant  of  the  surface-attached 
redox  species.  Laviron  theory  has  been  employed  by  others^’'*^  to  evaluate  electrode 
kinetics  from  cyclic  voltammograms  of  self-assembled  electroactive  monolayers.  We  have 
chosen  to  employ  a  Marcus  theory  formulation^^’^^  to  evaluate  rate  constants  from  AEpg^ 
results  since  Marcus  theory  is  more  appropriate  when  the  reaction  overpotential  is  not 
negligible  in  comparison  to  its  reorganization  energy.  A,.  The  pertinent  equation  is^*^ 
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Ked,ox  =  £ 


X  “*“71  ^ 
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l+exp(x) 


-dx 


(1) 


where  /i  is  the  electronic  coupling  (tunneling)  parameter,  p  the  density  of  electronic  states  in 
the  metal  electrode,  the  Boltzmann  constant,  T  the  Kelvin  temperature,  and  rj  the 
applied  overpotential.  This  equation  contains  an  integration  of  the  electronic  states  around 
the  Fermi  level;  we  refer^^^  to  it  as  Marcus-DOS  theory.  Equation  1  has  been  used^^  in 
digital  simulations  of  cyclic  voltammograms  to  examine  the  relationships  between  theoretical 
waveshapes  and  AEpg^  values  (AEp£^2  =  »?)»  and  reaction  and  experimental  parameters 
such  as  the  rate  constant  k“,  X,  T,  and  potential  sweep  rate,  v.  An  example  of  a  theoretical 
working  plot  of  AEpg^  vs.  log(v/k°),  calculated  for  125K,  is  given  in  Figure  2  as  is  an 
analogous  plot  based  on  Butler  Volmer  theory  (for  a  =  0.5). 

Figure  2  shows  that,  firstly,  AEp£^  values  predicted  by  Marcus  theory  are,  except 
when  AEpE^  is  relatively  small,  substantially  larger  than  those  predicted  by  standard  Butler 
Volmer  theory.  The  consequence  is  that  rate  constants  evaluated  when  AEpg^  is  large  with 
the  latter  electrode  kinetic  model  would  tend  to  be  too  small^"^®.  Secondly,  especially  for 
larger  X,  AEpE^K  values  are  relatively  insensitive  to  its  value  (although  waveshape  and  peak 
current  are  strongly  sensitive  to  A,,  vide  infra).  This  characteristic^’®  of  AEp^AK  insensitivity  to 
A,  is  beneficial  in  that  k°  can  be  estimated  reliably  from  AEpg^  without  accurate  prior 
knowledge  of  A,.  For  example,  in  the  present  work  an  experimental  cyclic  voltammogram 
for  a  CpFeCpC02(CH2)i2SH  monolayer  at  130K  and  50  mV/s  gave  a  AEpg^  =  410  mV; 
comparing  this  result  to  versions  of  Figure  2  calculated  for  130K  and  A.  =  0.6  v^.  1.0  eV 
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produces  k"  =  0.0009  vs.  0.00063  s‘\  respectively,  a  variation  of  only  ca.  50%  for  these  two 
assumptions  of  A,. 

AEpEAK  values  determined  as  a  function  of  potential  sweep  rate  at  different 
temperatures  for  CpFeCpC02(CH2)i2SH  monolayers  are  compared  to  theoretical  working 
curves  in  Figure  3.  The  theoretical  curves  were  calculated  for  X  =  0.7  eV.  The  figure  shows 
that  at  each  temperature  it  is  possible  to  select  a  rate  constant  that  produces  a  good  fit 
between  the  theoretical  curve  and  the  sweep  rate  data.  Table  I  gives  results  of  point-by¬ 
point  calculations  of  k“,  again  using  X  =  0.7  eV  working  curves.^®  Besides  the  temperature 
dependency  of  the  rate  constant  k°,  the  most  significant  aspect  of  Table  I  is  that  k“  does  not 
vary  with  v,  in  particular  k°  does  not  trend  downward  with  increasing  v.  A  potential  sweep 
rate  dependency  of  measured  rate  constants  is  generally  symptomatic  of  bias  of  AEpE^x 
values  caused  by  IRunc  problems;  its  absence  here  is  taken  as  evidence  that  the  resistance 
compensation  technique  employed  (see  Experimental)  satisfactorily  avoids  such  problems  in 
cyclic  voltammetry.  In  addition,  estimation  of  iRuNc  from  the  product  of  voltammetric  peak 
currents  and  directly  measured  Rung  (see  Table  I  footnote)  gives  voltages  which  are  small 
relative  to  the  measured  AEpg^  values,  even  before  positive  feedback  compensation  is 
applied.  Table  11  gives  k“  values,  averaged  over  different  sweep  rates,  for  a  larger  variety  of 
temperatures. 

The  variation  of  the  cyclic  voltammetrically  measured  k“  with  temperature  was 
analyzed  for  the  reorganization  energy  X  by  the  activation  relation 

k‘’  =  fipkj  exp  (-^) 


(2) 
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where  AG  the  free  energy  of  activation  is  X  /4  (assuming  on  the  basis  of  literature  data^  a 
negligible  entropy  of  activation).  Analysis  of  the  Table  n  k°cv  results  as  a  plot^^  of 

vs.  1/T  reveals  a  linear  activation  plot  (Figure  4(n))  with  slope  corresponding  to 
a  0.89  eV  reorganization  energy. 

The  result  of  Figure  4(0)  can  be  compared  to  the  dielectric  continuum  model^®  for 
barrier  control  by  the  energetics  of  solvent  dipole  (outer  sphere)  reorganization,  which  is 


X  = 


e'^N  11  11 

UL  (i—L)  (— --) 

Sire  2d'  ^e  e ' 

o  op  s 


(3) 


where  e  is  the  electronic  charge,  N  Avogadro’s  number,  the  permittivity  of  free  space,  a 
the  reactant  radius  and  d  its  separation  from  the  electrode  surface,  and  e^p  and  the 
optical  and  static  dielectric  constants,  respectively.  Several  assumptions  are  involved  here: 

(/)  Reactions  of  ferrocene  species  are  generally  thought  to  be  mainly  outer  sphere  barrier- 
controlled^,  which  is  assumed  here,  (zi)  d  is  18.7x10'^°  m  based  on  an  all  trans  30°  tilted 
alkane  chain  configuration‘s  and  treating  the  ester  hnkages  as  equivalent  to  two  additional 
methylene  units.^^  (Hi)  The  small,  linear  temperature  dependencies^’^  of  e^p  and  for  EtCl 
and  PrCN  have  been  extrapolated  to  our  range  of  temperatures^  and  combined  as  an 
"average"  dielectric  continuum  with  "average"  dielectric  properties  e^p  and  that  are  sums 
of  the  mole-fi'action-weighted  individual  constants.  The  sensitivity  of  X  to  temperature  is 
very  small  (predicted  X  is  0.78  eV  at  293K  and  0.76  eV  at  175K),  and  since  the  more 
important  e^p  term  is  similar  for  EtCl  and  PrCN,  an  averaged  dielectric  seems  reasonable. 
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(rv)  The  ferrocene  is  actually  positioned  at  a  high  polarity  (EtCl/PrCN)-low  polarity  (the 
alkane  chains)  interface,  whereas  our  use  of  Equation  3  assumes  uniform  solvation  by  the 
high  polarity  medium.  An  analysis^  of  solvent  spatial  correlation  effects  may  be  a  preferred 
approach.  That  Equation  3  is  a  simplification  of  the  actual  situation  has  also  been  noted  by 
previous  workers^. 

Irrespective  of  the  assumptions  noted  above,  the  reorganization  energy  X  =  0.89  e V 
derived  from  the  activation  plot  in  Figure  4(0)  is  only  slightly  larger  than  the  0.76  eV 
prediction  of  Equation  3.  For  reasons  described  below,  there  are  difficulties  in  other 
approaches  to  measuring  X,  namely  voltammetric  waveshape  analysis  and  potential  step 
experiments,  that  may  produce  less  reliable  results  for  ionically  resistive  solutions  and 
kinetically  disperse  monolayers. 

Extrapolation  of  the  Figure  4(n)  activation  plot  gives  a  predicted  room  temperature 
^273°  =  36  s'^  and  a  product  /ip  =  6.5x10^  eV'^s'^  The  room  temperature  rate  constant  is  not 
actually  measureable  in  the  EtCl/PrCN  solvent  since  the  self-assembled  monolayers  are 
stable  in  it  only  at  low  temperatures.  The  data  can,  however,  be  compared  to  room- 
temperature  k°  results  in  aqueous  acid  by  Chidsey^’^'’  of  ca.  40  s'^  and  2.1x10^  eV  ^s  ^  for  k° 
and  MP»  respectively.  The  agreement  between  the  k"  and  /xp  values,  considering  the 
difference  in  medium  and  the  kinetic  dispersion  {vide  infra)  of  the  present  monolayers,  is 
quite  reasonable. 

As  noted  above,  we  have  shown^"^®  that  cyclic  voltammetric  waveshapes  are  in 
Marcus-DOS  theory  quite  sensitive  to  the  value  of  reorganization  energy  X.  Waveshapes 
become  broadened  when  X  is  not  much  larger  than  AEpg^,j^2.  The  waveshape  sensitivity 
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can,  however,  be  employed  to  measure  X  onfy  when  the  experimental  voltammetry  does  not 
also  experience  broadening  from  other  sources,  such  as  a  distribution  of  formal  potentials 
(E°’)  or  of  the  actual  rate  constants  among  the  surface  redox  site  population.  We  refer  to 
these  surface  population  inhomogeneity  effects,  collectively,  as  kinetic  dispersion^.  The 
presence  of  kinetic  dispersion  wiU  cause  a  X  value  derived  from  a  Marcus-DOS  waveshape 
analysis  to  be  erroneously  small,  as  was  demonstrated  in  a  recent  comparison^’^®  of 
voltammetry  from  one  CpFeCpC02(CH2)i6SH  monolayer  which  did  not  exhibit  kinetic 
dispersion  (as  established  by  linear  ln[i]  vs.  time  plots  in  potential  step  experiments^)  with  a 
different  monolayer  sample  which  by  the  same  criterion  was  kinetically  disperse.  This 
comparison^’^®  also  gave  the  important  result  that  the  rate  constants  k°  obtained  from  the 
kinetically  homogeneous  and  kinetically  disperse  monolayers  by  use  of  working  plots  such  as 
Figure  2  were  in  good  agreement.  This  agreement  suggests  that  the  averaging  involved  in 
determining  a  rate  constant  with  Figure  2  does  not  significantly  bias  the  result  towards  the 
faster  or  the  slower  reacting  sites  among  the  ferrocene  surface  population.  Calculations^®'’ 
are  also  consistent  with  the  experimental  observation^^®,  indicating  only  a  two-fold,  nearly 
temperature  independent  upward  bias  in  a  cyclic  voltammetric  k°  measurement  of  a 
kinetically  disperse  monolayer. 

AU  ferrocene  alkanethiol  monolayers  that  we  have  examined^'”’^^’^®’^^®  in  the  low 
temperature  EtCl/PrCN  solvent  exhibit  substantial  broadening  of  cyclic  voltammetric 
waveshapes  and  (vide  infra)  non-linear  ln[i]  vs.  time  plots  in  potential  step  experiments.  The 
kinetic  dispersity  of  these  monolayers  occurs  probably  in  part  by  intercalation  of  the  organic 
solvent  into  the  monolayer.  Figure  5  shows  a  comparison  of  waveshapes  for  voltammetry  of 
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CpFeCpC02(CH2)i2SH  monolayers  at  different  temperatures,  to  waveshape  calculations 
based  on  Marcus-DOS  theory  (Equation  1).  It  is  clearly  possible  to  select  theoretical 
parameters  that  produce  reasonable  fits  to  the  experimental  results,  and  the  rate  constants 
obtained,  k“  =  0.0018,  0.014,  0.09  s'^  at  130,  150,  and  170K,  are  not  very  different  from 
those  obtained  using  Figure  3  with  the  AEp^AK  values  of  these  voltammograms,  0.0008,  0.011, 
and  0.088  s'^,  respectively.  However,  the  reorganization  energies,  X  =  0.29,  0.29,  and  0.18 
eV  at  130,  150,  and  170K,  respectively,  that  are  required  to  fit  the  theoretical  waveshapes  to 
the  experimental  ones  are  very  low.  These  latter  reorganization  energy  results  are  beUeved 
to  be  unreliable  and  are  presented  solely  to  illustrate  a  pitfall  of  kinetic  dispersion. 
Measurement  Using  Potential  Steps.  Chidsey^  and  Finklea^  have  measured  heterogeneous 
electron  transfer  rate  constants  in  electroactive  self-assembled  monolayers  using  current¬ 
time  transients  produced  by  potential  steps  to  various  overpotentials.  Electrolysis  of  the 
monolayer  is,  ideally,  a  first  order  process,  giving  a  current  transient  described  by  the 
relation 

i  =  k  Q  exp(-k  f)  (4) 

T}^  app^  ■*  '  app,V  ' 

where  the  sum  of  the  forward  and  reverse  electron  transfer  rate  constants  at 
overpotential  rj  (E^pplied  ‘  E*”)  is  given  by  the  slope  of  a  ln[i]  vs.  time  plot,  and  Q  is  the 
total  charge  passed  in  the  electron  transfer  reaction.  The  linearity  of  such  a  plot  of 
Equation  4  is  an  indicator  that  the  reacting  surface  population  is  kinetically 
homogeneous.^*^ 

Figure  6  shows  ln[i]  vs.  time  plots  acquired  at  140K  by  stepping  from  E°’  by  +300, 
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+ 150,  and  +70  mV,  using  positive  feedback  iR  compensation^.  Immediately  obvious  are 
the  substantial  differences  in  current  and  time-scale  as  f]  is  varied,  and  that  the  plots  exhibit 
curvature  toward  smaller  slope  at  all  overpotentials.  The  latter  effect  is  symptomatic  of  a 
kinetically  disperse  or  non-uniform  monolayer.  We  have  analyzed  these  dispersive  kinetics 
by  dividing  the  current-time  transient  into  segments  of  charge  passed,  which  serves  to  define 
the  different,  sequentially  reacting,  populations  of  surface  ferrocene  sites.  The  vertical  Hnes 
on  Figure  6  indicate,  for  example,  the  beginning  and  end  of  passage  of  charges  which  are 
the  0.5  to  0.8  fraction  of  the  total  reaction  charge  (nFA[r]/2).  The  plots’  slopes  during 
passage  of  each  such  increment  of  charge  (or  of  any  smaller  increment)  represents  the 
kinetics  of  that  particular  segment  of  the  reacting  surface  population.  The  rate  constants  of 
such  population  segments  can  then  be  compared  as  to  their  dependency  on  the  applied 
overpotential. 

Figure  7  plots  rate  constants  k^pp,,,  for  the  successive  0.5  to  0.6,  0.6  to  0.7,  and  0.7  to 
0.8  charge-based  segments  of  reacting  population  in  experiments  in  which  the  size  of  the 
potential  step  (rj)  was  varied,  and  for  experiments  at  several  temperatures.  These  results 
show  that  the  variation  of  rate  constants  for  the  three  charge  segments  (the  three  over¬ 
lapping  data  points  at  each  rj  and  temperature)  is,  while  noticeable,  small  in  comparison  of 
the  variations  of  rate  constant  due  to  changes  in  rj  or  in  temperature.  We  see  that  the  rate 
constants  in  Figure  7  do  not  follow  the  exponential  increase  with  overpotential  classically 
expected  from  Butler  Volmer  theory,  but  instead  fold  over  at  higher  rj  as  anticipated  by 
Marcus-DOS  theory^’^  when  rj  is  not  negligible  with  respect  to  reorganization  energy  X.  It 
is  further  significant  that,  like  results^  for  kinetically  uniform  ferrocene  monolayers  in 
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aqueous  acid,  the  variation  in  rate  constant  with  rj  is  symmetrical  about  r]  =  0.  As 
discussed  by  Finklea^  such  symmetry  is  evidence  for  a  through-bond  tunneling  mechanism. 

Over-lying  each  of  the  sets  of  log[k^P,,N]  vs.  rj  data  in  Figure  7  are  best  fit 

cui-ves( - )  derived  by  comparison  of  Equation  1  to  the  data.  Table  II  contains  the  values 

of  k°  and  X  corresponding  to  the  fits  in  Figure  7  and  for  further  experiments.  Examining 
these  results,  we  see  firstly  that  rate  constants  derived  from  the  potential  step  analysis  are  in 
generally  good  agreement  (within  a  factor  of  two)  with  those  obtained  (vide  supra)  from 
cyclic  voltammetry.  Secondly,  the  reorganization  energy  obtained  from  the  log[kAPP,,7]  V5.  rj 
fits  of  Equation  1  at  the  highest  temperature,  0.80  eV  at  150K,  is  close  to  that  obtained,  0.89 
eV,  from  the  activation  analysis  of  the  k’cv  data  (Figure  4(0)).  Unexpected,  however,  is  the 
apparent  temperature  dependence  of  X  of  the  reorganization  energy,  in  which  we  see  X 
decrease  to  0.46  eV  at  the  lowest  temperature.  Thirdly,  the  pre-factor  jup  of  Equation  1 
obtained  in  the  fitting  is  at  150K,  5.3x10^  eV^s'\  This  pp  value  can  be  compared  to  data  by 
Chidsey^  which  gives  pp  =  2.1x10®  eV'^s'^  for  CpFeCpC02(CH2)i2SH  monolayers  in  aqueous 
HCIO4  at  273K  Fourthly,  below  150K,  consistent  with  the  changes  seen  in  X,  the  values  of 
pp  (not  shown)  corresponding  to  the  Table  II  data  become  much  smaller,  falling  to  4.2x10^ 
eV^s’^  at  125K  As  explained  next,  we  believe  the  X  results  at  the  lower  temperatures 
(shown  parenthetically  in  Table  H)  are  not  reliable. 

The  rate  constants  for  the  0.5  to  0.8  population  segments  were  selected  for  analysis  in 
Figure  7  since  these  data  are  potentially  less  distorted  by  residual  uncompensated 
resistances  than  those  collected  during  earlier  reacting  segments  where  the  currents  are 
larger.  One  consequence  of  a  resistance  effect,  if  present,  would  be  to  depress  k^p  values 


/ 


16 

obtained  at  larger  rj  (where  higher  currents  flow)  more  than  those  obtained  at  lower  ry;  the 
added  curvature  thus  induced  in  log[k^p]  vj.  rj  plots  would  consequently  depress  the  result 
for  X.  This  error  would  be  most  serious  at  lower  temperatures  since  ionic  conductivities 
decrease  with  decreasing  temperature.  (We  do  not  on  the  other  hand  believe  that 
resistance  effects  are  a  problem  in  the  cyclic  voltammetric  experiments  due  to  both  the 
results  of  Table  I  and  because  the  currents  passed  in  those  experiments  are  much  smaller 
than  in  Figure  6.) 

Even  were  resistance  effects  absent  in  the  potential  step  experiments,  kinetic 
dispersity  in  the  monolayers  can  cause  spurious  results  in  data  analyses  Uke  Figure  7.  As 
described  elsewhere^*’,  simulation  of  the  effects  of  a  spread  in  ferrocene  formal  potentials 
E°’  on  a  potential  step  rate  measurement,  shows  that  the  E°’  dispersion  causes  (/)  curved 
ln[i]  vs.  time  plots,  and  (ii)  depressed  values  for  k",  X,  and  jup  upon  analysis  with  Equation  1, 
for  all  charge  segments  and  especially  early  ones.  We  are  as  a  consequence  skeptical  of  the 
parenthetically  shown  potential  step  results  in  Table  n  for  A,. 

The  k*  results  from  potential  step  experiments  are  plotted  according  to  Equation  3  in 
Figure  4(0),  and  give  from  the  slope  X  =  0.66  eV.  This  value  of  the  reorganization  energy 
is  smaller  than  that  obtained,  0.89  eV,  from  the  cyclic  voltammetric  data.  Figure  4(0)  and 
slightly  smaller  than  predicted  by  theory,  0.76  eV.  The  k"  values  determined  in  potential 
step  experiments  (and  the  X  values  derived  from  them)  are  apparently  somewhat  less 
affected  by  the  problem  noted  above. 

In  summary,  the  resistive  low  temperature  environment  and  the  kinetic  dispersion  in 
the  CpFeCpC02(CH2)i2SH  monolayers  present  a  challenging  measurement  situation  within 
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which  are  pitfalls  that  we  believe  can  lead  to  distortions  in  the  kinetic  results.  The 
derivation  of  rate  constants  k”  from  analysis  of  cyclic  voltammetric  AEpg^  results  and  of  A 
by  activation  plots  of  those  rate  constants  appear  to  be  most  immune  to  such  distortions. 
ACKNOWLEDGEMENT.  This  research  was  ftmded  in  part  by  grants  from  the  Office  of 
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temperature  voltaimnetric  waves.  Further,  calculations  for  potential  step  experiments 
show^  that  a  spread  among  the  thermodynamic  potentials  can,  since  applied 
overpotential  r?  is  referenced  to  the  average  E°’,  cause  curved  ln[i]  vs.  time  plots  for 
such  experiments,  thus  mimicking  a  dispersion  of  rate  constants  among  the  surface 
redox  sites.  This  apparent  dispersion  would  add  to  any  real  dispersion  among  rate 
constants  that  might  also  be  present.  The  E°’ -based  kinetic  dispersion  also  has  a 
substantial  effect  on  plots  of  log[k^p]  vx.  rj,  introducing  curvature  at  large  rj  that  in 
analysis  of  such  plots  by  Equation  1  causes  smaller  apparent  values  of  k°,  X,  and  up. 
(b)  G.  K.  Rowe,  M.  T,  Carter,  J.  N.  Richardson,  R.  W.  Murray,  submitted  for 
publication. 

27.  The  transients  of  Figure  6  were  corrected  for  background  currents  with  transients 

resulting  from  with  the  same  magntiude  potential  step  in  the  double  layer  region  that 
lies  negative  of  the  surface  CpjFe"^/®  couple. 


TABLE  I.  Heterogeneous  electron  transfer  rate  constants  obtained  from  cyclic  voltammetry 
as  a  function  of  potential  sweep  rate  and  temperature.® 


a)  Uncompensated  resistances  determined  using  AC  impedance  at  T  =  125,  130,  135,  145, 
and  150K  are  Rung  =  5.5x10^,  2.3x10^,  9.3x10^,  2.7x10^,  and  1.6x10^  n,  respectively.  These 
are  values  before  positive  feedback  compensation  is  applied,  b)  from  comparison  of 
experimental  AEpg^  to  those  of  digitally  simulated  voltammograms,  assuming  X  =  0.7  eV 
(see  text). 


23 


TABLE  n.  Rate  parameters  from  analysis  of  log(k3pp)  vs.  rj  plots  and  cyclic  voltammetry  as 
a  function  of  temperature. 


T(K) 

k°,  s'^ 

(E  step)® 

X,  eV 
(E  step)® 

k°,  s-^ 

(CV)” 

125 

0.0005 

(0.46) 

0.00029 

130 

0.001 

(0.48) 

0.0008 

135 

0.002 

(0.49) 

0.0017 

140 

0.003 

(0.56) 

0.0032 

145 

0.004 

0.78 

0.0064 

150 

0.008 

0.80 

0.012 

155 

— 

— 

0.020 

160 

— 

— 

0.034 

165 

— 

— 

0.057 

170 

— 

— 

0.090 

175 

— 

— 

0.136 

a)  From  analysis  of  log(k^p  vs.  rj  plots  (see  text).  Activation  analysis  of  these  data  gives  X 
=  0.66  eV,  MP  =  7.4x10^  eV’^s’\  and  an  extrapolated  k273  =  3.2  s'\  all  smaller  than  data 
from  CV  analysis,  b)  from  comparison  of  experimental  AEpg^  to  those  of  digitally 
simulated  voltammograms,  assuming  X  =  0.7  eV  (see  text).  Activation  analysis  of  these  data 
gives  X  =  0.89  eV  and  mP  =  6.5x10^  eV'V^ 
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FIGURE  CAPTIONS 

Figure  1.  Cyclic  voltammograms  (50  mV/s)  of  mixed  CpFeCpC02(CH2)i2SH/ 

CH3(CH2)nSH  monolayer  in  2:1  (v:v)  EtO/PrCN  solvent  at  T  =  130K  ( - ),  150K  (--- 

-),  and  170K  ( — ).  Electroactive  ferrocene  coverages  are  1.4x10'^°,  1.5x10'^°  and  1.6x10"^° 
mol/cm^,  respectively. 

Figure  2.  AEpe^k  log(v/k°,  units  are  volts)  at  125K  from  digitally  simulated  cyclic 
voltammograms  based  on  Marcus-DOS  theory  (Eqn.  1)  as  a  function  of  reorganization 

energy:  X  =  0.2  eV  (□),  0.4  eV  (V),  0.7  eV  (A),  and  1.0  eV  (  +  ),  ( - ^)  represents  result  using 

Butler-Vohner  theory  for  a  =  0.5. 

Figure  3.  Experimental  AEpE^K log(v/k“  units  are  volts)  at  130K(n)  and  k°  =  0.0008  s\ 
140K(O)  and  k"  =  0.0032  s‘\  and  150K(A)  and  k*  =  0.012  s'^,  compared  to  digitally 
simulated  cyclic  voltammograms  at  those  temperatures,  based  on  Marcus-DOS  theory 
(Equation  1)  and  a  reorganization  energy  A,  =  0.7  eV. 

Figure  4.  Activation  plot  of  Equation  2  for  experimental  k°  data  in  Table  n,  from  cyclic 
voltammograms  (□)  and  log[k^p^]  vs.  r\  plots  (O),  slopes  and  intercepts  give  X  =  0.89  eV 
and  0.66  eV,  and  mP  =  6.5x10^  eV'^s'\  and  7.4x10^  eV'^s*\  respectively,  were  calculated  from 
the  slopes  of  these  plots. 

Figure  5.  Waveshape  fits  of  experimental  (oxidation  branch)  voltammetric  waves,  v  =  50 
mV/s  at  130K(n),  150(0),  and  170K(A),  to  waveshapes  calculated  based  on  Equation  1  and 
k“  =  0.0018  s-\  X  =  0.29  eV  (130K);  k”  =  0.014  X  =  0.29  eV  (150K);  and  k“  =  0.09  s'^  X 
=  0.18  eV  (170K). 

Figure  6  ln[i]  vs.  time  transients  obtained  at  140K  at  rj  =  +300  mV  (panel  A),  +150  mV 
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(panel  B),  and  +70  mV  (panel  C).  Vertical  lines  represent  boundaries  of  the  fractions  0.5 
to  0.8  of  the  total  charge  passed  in  the  reaction,  nFAr/2.  The  average  kgpp  ^  over  this 
charge  segment  is  30,  0.63,  and  0.036  s'^  for  the  300,  150,  and  70  mV  potential  steps, 
respectively. 

Figure  7.  (O)  Oxidation  and  reduction  branches  of  a  WOK  composite  log(k^p ,,)  vs.  rj  plot 
(Tafel  plot)  constructed  using  values  from  0.5-0.6,  0.6-0.7  and  0.7-0.8  charge  segments 
(bottom-to-top  of  three  O  shown)  of  total  reacting  charge  in  the  step;  (□)  same  but  130K 

and  oxidation  branch  only;  (A)  same  but  150K  and  oxidation  branch  only;  ( - )  are  best  fit 

lines  from  Equation  1  for  parameters  (130K)  k“  =  0.001  s‘\  k  =  0.48  eV,  and  /ip  =  8.2x10^ 
eV-^s-\  (WOK)  k“  =  0.0025s*\  k  =  0.52  eV,  and  mP  =  2.0x10^  eV-^s-\  and  (150K)  k“  =  0.008 
s-\  X  =  0.80  eV,  and  MP  =  5.3x10^  eV^'^ 
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